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Abstract 

We have considered a numerical scheme for the calculation of the equilibrium properties 
of spin-^ XY chains. Within its frames it is necessary to solve in the last resort only the 
2N X 2N eigenvalue and eigenvector problem but not the 2^ x 2^ one as for an arbitrary 
system consisting of N spins ^. To illustrate the approach we have presented some 
new results. Namely, the xx dynamic structure factor for the Ising model in transverse 
field, the density of states for the isotropic chain with random intersite couplings and 
transverse fields that linearly depend on the surrounding couplings, and the zz dynamic 
structure factor for the Ising model in random transverse field. The results obtained are 
hopped to be useful for an interpretation of observable data for one-dimensional spin-^ 
XY substances. 
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Statistical mechanics of the one-dimensional spin-| XY models has a history going back 
more than three decades to the pioneering paper by Lieb, Schultz and Mattis [1] in which 
it was noted that such systems are as a matter of fact the non-interacting spinless fermions 
and therefore, a lot of statistical mechanics calculations can be performed exactly. Although 
many papers concerning thermodynamics, spin correlations and their dynamics have appeared 
since that time, several problems still call for efforts. One of them regards the study of spin 
dynamics: in contrast to the time- dependent correlation function between ^-components of 
two spins [2], the xx time-dependent correlation function has been derived exactly only in 
some limiting cases (T = 0, the critical value of transverse field, T = oo) [3-9]. Some recent 
attempts to calculate this correlation function arc presented in Rcfs. [10-14]. New difficulties 
arose when the one-dimensional spin-| XY models were started to be discussed in random 
versions. The analytical results obtained here are not so impressive as for the perfect case and 
they are restricted to special types of disorder [15-20] or renormalizat ion-group analysis [21]. 
A discovery of the quasi-one- dimensional spin-| systems (Cs(Hi_xE)x)2P04, PbHi_xDxP04, 
PrCls, CsCuCls, CsCui-xMnxCls, J-aggregates etc.) gave arise to additional interest in such 
calculations since a lot of data obtained in the dynamic experiments for such materials are 
available. Some approximate studies of dynamic properties were inspired by corresponding 
measurements [22-27], however, such estimations contain uncontrolled mistakes. Additional 
interest in such calculations was caused by the very recent studies of quantum phase transitions 
in disordered systems [21,28-30]. 

Since there are notorious difficulties in the analytical study of some properties of the spin-i 
XY chains it is naturally to try to obtain the desired results numerically. Few earlier attempts 
[31,32] faced the 2^ x 2-^ eigenvalue and eigenvector problem that restricted the computations 
to N 10. However, a peculiarity of spin-| XY chains provides an evident possibility to 
compute all equilibrium quantities facing in the last resort only the 2N x 2N eigenvalue and 
eigenvector problem [33-35]. This fact allows to consider rather long chains {N ~ 100 — 10000) 
and therefore to study reliably the dynamics of spin correlations or the influence of disorder 
on observable properties. In this paper we shall briefly explain such a numerical approach for 
calculation of the equilibrium properties of the spin-i XY chains. Besides we shall demonstrate 
how does the approach work computing for this purpose the xx dynamic structure factor of 
the transverse Ising model and some thermodynamic and dynamic properties of several random 
XY models. 

We shall consider N spins one-half arranged in a row with the following Hamiltonian 

N N-1 
j=l j=l a,P=x,y 

where VLj is the transverse field at site j and Jj^ is the interaction between a and j3 spin 
components at the sites j and j + 1. Introducing instead of the spin raising and lowering 
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operators via the Jordan- Wigner transformation the Fermi operators Cj, Cj one finds that 
the Hamiltonian (1) is a bihnear fermion form that can be put into diagonal form by hnear 
transformation r]'^ — J^jLi (j^lj'^j + dlj'^t) ■ Similarly to [1] it can be shown that if 



(gfc, hfe) 




Afc(gfc,hfe) 



(2) 



with 



Jjji—ly 



Bij = 



Evidently, the knowledge of 



Jj^)] the Hamiltonian transforms into H 



AfcS or their distribution p{E) = Z^^i S{E—Ak) yields thermodynamics of the spin system (1). 
For the calculation of spin correlation functions it is convenient to introduce the auxiliary oper- 
ators ip+ = ct + Cj = EjLi (%jVp + %jVi^ and ifj = c+ - Cj = EJLi {^pjVp - %jVp), where 
%j = 9pj+hpj^ = 9pj-hpj- Since s] = ^(pt^I ■ ■ ■ (ft-iy^J-ift^ 4 = ^^'/'i^^F • • ■ ^t-i^fij-i^pj 
and s5 



.J the calculation of thermodynamic average of a product of spin operators 

reduces to exploiting the Wick-Bloch-de Dominicis theorem with the result that is in fact the 
Pfaffian of corresponding antisymmetric matrix constructed from the elementary contractions 



N 



p=l 



N 



p=l L 



^pj ^pm 



^(Ap) ^(-Ap,j 



{^j{t)0 = E 

p=l 

N 

-E 



\|/* dj 



:F{Ap) T{-Kp) 



with ^(x) = (^1 + e'^'^j e~^^*. Thus, the solution of the 2N x 2N eigenvalue and eigenvector 
problem (2) completely determines thermodynamics, spin correlations and their dynamics for 
the model (1). 

In what follows we have collected several new results obtained within the frames of the 
described approach. We shall start from the xx dynamic structure factor Sxx{i^^ ^) = Yln=i 6^**" 
dte~^'*'e^'^*(s^(t)sj^„) for the uniform transverse Ising chain (Qj = Q = 0.2, Jj^ = J = — 1, 
Jj^ = Jj^ = Jj^ = 0). At first we solved the eigenvalue and eigenvector problem (2) for 
N — 280 obtaining in result A^, ^kj, ^^j. Then for several values of temperature we calculated 
the required elementary contractions and computed the relevant Pfaffian 
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At last we performed the integration over the time up to t = 1600 with e — 0.001 and took into 
account in the sum over n up to 100 neighbours. The obtained dependences Sxxii^, vs. cu for 
different k at low temperature and SxxiOjOu) vs. to for few temperatures are depicted in Fig. 1. 
As it can be seen from Fig. 1 Sxx{i^, exhibits two peaks. At (3 — b and k — Q one finds a high 
peak at a;i = 0.00 and a low and broad one at a;2 = 0.76. As k increases the height of the first 
peak decreases and it shifts towards high frequencies {uji ^ 0.11 at k = |, 0.25 at k = |, 
LOi 0.34 at K = oil 0.38 at k = tt), whereas the width of the second peak decreases, 
its height increases and it moves towards high frequencies {u>2 ~ 0.80 at k — uj2 ^ 0.88 at 
K — ^, UJ2 ^ 0.96 at re = ti;2 ~ 1-01 at K = tt). As the temperature increases Sxx{0,'^) 
qualitatively remains the same: the heights of both peaks decrease and the high-frequency peak 
shifts slightly towards high frequencies {uj2 ~ 0.76 at /9 = 5, c<j2 ~ 0.77 at P = 4, 002 ^ 0.81 
at f3 = 3, 002 ^ 0.86 at /3 = 2, ~ 0.91 at fj = 1, uj2 ^ 0.95 at fi = 0.1, 002 ^ 0.96 at 
f3 = 0.001). The discussed case of small transverse field may be of interest for understanding 
the dielectric measurement and neutron scattering data for quasi-onc-dimensional hydrogen- 
bonded ferroelectrics. However, detailed comparison with experimental results demands an 
introducing of weak interchain interactions (see [24,26]) that requires a separate study. 

Next example regards to thermodynamics of the isotropic XY chain {Jj^ = Jj^ = Jj, 
Jj^ = Jj^ = 0) with the random intersite couplings given, for example, by the Lorentzian 
probability distribution density p{. . . , Jj , . . .) = IljLi - jj-^^^jr^ and the transverse fields that 
depend linearly on surrounding couphngs f2j — — |(Jj-i + Jj — 2Jo), | a |> 1. In Ref. [20] 
the random- averaged density of states p{E) was calculated that yielded thermodynamics of the 
model. In particular, it was shown that the introduced randomness may cause the appearance 
of non-zero averaged transverse magnetization rn^ — Jf^ dEp{E) tanh ^ at zero averaged 
transverse field = {p{E) is not symmetric with respect to the change E — Qq ^ —{E — Qq) 
and at r = -2m; = dEp{E) - dEp{E) y^O). We considered the described model 
for A'" = 15000 with Jq = 1 for a = ±1.01 with the Lorentzian and Gaussian random couplings 
(the probability distribution for the latter case reads p{. . . , Jj, . . .) — Ilf=i -;j=^ cxp —^-^^^^^r- ) 
and presented the results of computations for one random realization in Fig. 2. The displayed 
results are in complete agreement with the ones derived analytically for the Lorentzian disorder. 
However, the developed numerical procedure allows to study an arbitrary disorder. The per- 
formed calculations have indicated the appearance of 7^ for = at low temperature for 
the correlated disorder. This can be clearly demonstrated by the results of calculation of a num- 
ber of negative and positive eigenvalues A^ (2) M- and since at T = —2rfT^ = ■ 
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Putting a — 1.01 for a certain Gaussian random realization with a — 0.25(1) that yields 
wEf=i Jj = 0.999757 (0.999027) we found jV- = 7192(6024) and Af+ = 7808(8976). Another 
Gaussian random reahzation with a = 0.25(1) that yields ^ EjLi Jj = 1-000118 (1.000473) gave 
jV- = 7187(6073) and Af+ = 7813(8927). These results definitely point out the appearance of 
"spontaneous magnetization" due to disorder. 

We end up with the calculation of the zz dynamic structure factor S'^^(k, cu) = Yln=i ^^'^^ I^oo 
dte-^l^lei"^* ((s|(t)s|+„) - (s|)(s|+„)) for the Ising chain (Jf^ = J = 1, Jf^ = Jf = Jf = 0) 
in the random transverse field defined by the probability distribution density p{. . . ,VLj, . . .) = 
YljLi [x6{^lj) + (1 — x)5{ilj — 0.5)], < X < 1. We computed the correlation functions 

"^lOO+n n ) - (V'l00(0V^100+n)(<^100(i)<^100+n) + {A{t) 

{fiooi'^)vtoo+n) for 250 random chains of 200 spins, performed the integration over time t with 
e = 0.005 and the summation over neighbours n. The obtained random-averaged zz dynamic 
structure factor that for = is presented in Fig. 3. The depicted plots demonstrate how the 
frequency-dependent zz structure factor rebuilds from the Ising type behaviour to the trans- 
verse Ising type behaviour as the concentration of sites with transverse field increases from 
(x = 1) to 1 (x = 0). The obtained dynamic structure factor exhibits a lot of structure that 
is induced by the disorder arrangement of two values of transverse field and 0.5. It appears 
that each well-defined peak for small concentrations of 1 — a; is connected with ^'^^(O, a;) for a 
certain chain determined by local environment of spin at j — 100 (for example, . . . 000 . . . , 
...OQOO..., ...0120..., ...OnnOO..., ...OnnO..., ...Onon... etc., the transverse field 
Qioo is written in bold font, not written Q^s do not infiuence S^ziOjCo)). With decreasing of x 
a number of possible local structures in the vicinity of j — 100 (and thus a number of peaks) 
increases and the peaks appear almost at all frequencies. However, the difference in their 
heights is conditioned by the probability of their appearance that is large. As a result one 
gets fine structure that transforms into the smooth curve only in the limiting case x — 0. The 
described random model has a simple interpretation in connection with partially deuterated 
quasi-one-dimensional hydrogen-bonded ferroelectrics. However, a study of relevant there xx 
spin dynamics is more cumbersome and will be reported separately. 

In summary, we have presented the numerical approach suitable for calculation of time- 
dependent correlation functions for non-random and random spin-| XY chains. We have 
illustrated the numerical procedure deriving some new results. It is relevant to mention here 
the papers [36,37] devoted to numerical calculations of the xx time-dependent spin correlation 
functions for the isotropic XY model. The authors used explicit expressions for the elemen- 
tary contractions (in contrast to our formulae that also fit the random models) and computed 
determinants of corresponding antisymmetric matrices that yielded only the square of corre- 
lation functions (this causes some difficulties in further calculation of the dynamic structure 
factor or susceptibility as well as in the study of random models). We hope that the numerical 
analysis of the properties of spin-| XY chains will be useful for understanding the results of 
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the corresponding measurements on quasi-one-dimensional hydrogen-bonded ferroelectrics hke 
CSH2PO4, Cs(Hi_xDx)2P04, PbHP04, PbHi_xDxPC)4 (neutron scattering, dielectric mea- 
surement) and J-aggregates (absorption and emission spectra). However, a comparison of 
theoretical predictions and experimental data requires further studies. 
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List of figure captions 



FIG. 1. The dynamic structure factor for the Ising chain in transverse field, (a) Sxxi,!^^^) vs. 
uj (3 — b lor K — 0, |,|,^,7r (curves 1,...,5, respectively) and (b) S^xi^,^) vs. u for 
/? = 5, 4, 3, 2, 1, 0.1, 0.001 (curves 1, . . . , 7, respectively). 



FIG. 2. The density of states and transverse magnetization for the isotropic XY chain with 
random couplings and transverse fields that depend linearly on surrounding couplings. p{E) 
vs. E — VLq (a,c) and — vs. VLq aX [5 = 100 (b,d) for the Lorcntzian (short dashed lines) and 
Gaussian (solid lines) disorder with V = a = 0.25 (a,b) and F = a = 1 (c,d); the results for 
non-random case F = cr = are depicted by long dashed lines. 



FIG. 3. The frequency dependence of the averaged dynamic structure factor 5*22(0,0;) for the 
Ising chain in random transverse field at P — 5. 
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